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Abstract. The large body of near infrared observations 



presented in Testi et al. (1997; 1998) are analysed with 
the aim of characterizing the young stellar clusters sur- 
rounding Herbig Ae/Be stars. The results confirm the ten- 
dency of early Be stars to be surrounded by dense clus- 
ters of lower mass "companions" , while Ae stars are never 
found to be associated with conspicuous groups. The tran- 
sition between the different environments appears to occur 
smoothly from Ae to Be stars without a sharp threshold. 

No correlation of the richness of the stellar groups de- 
tected is found with the galactic position or the age of the 
central Herbig Ae/Be star. 

The stellar volume densities estimated for the groups 
surrounding pre-main-sequence stars of intermediate mass 
show the transition from the low density aggregates of T 
Tauri stars and the dense clusters around massive stars. 

Only the most massive stars (10-20 M©) are found 
to be associated with dense (~ 10 3 pc -3 ) stellar clusters. 
This is exactly the mass regime at which the conventional 
accretion scenario for isolated star formation faces theo- 
retical problems. Thus our findings strongly supports the 
idea that the formation of high-mass stars is influenced by 
dynamical interaction in a young cluster environment. 

Key words: stars: formation - stars: pre-main-sequence 
- infrared: stars - open clusters and associations: general 



1. Introduction 

The tendency of stars to form in approximately coeval 
groups, rather than in isolation, is b y now a well estab- 
lis hed o bservational f act (L ada et al. 1993[ Zinnecker et 



al. |l993t G omez et al. |l993| ; Hillenbrand |l995i |l997fc Testi 



et al. 1997 ). A result which sets stringent constraints to 
the theory is the fact that the "richness" of the groups de- 
pends on the mass of the most massive star in the cluster: 
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low-mass stars form in loose groups of a few objects per 
cubic parsec (Gomez et al. 1993), which in the following 



we will call "aggregates" , while high-mass stars are usu- 
ally found to form in dense clusters with densities up to 
10 4 objects per cubic parsec in the case of t he Tr apezium 
cluster (see e.g. McCaughrean & Stauffer 1994 ; Hillen- 
brand 1997 ). The transition between these two modes of 
formation occurs in the mass interval 2 < M/M© < 15. 



Herbig Ae/Be stars (Herbig 1960 ) are pre-main-sequence 
(PMS) stars of intermediate-mass located outside complex 
star forming regions. These stars are sufficiently old (age 
~ 0.5 — 5 Myr) to be optically visible, but still young 
enough that any population of lower mass stars born in the 
same environment has not had time to move away from 
their birthplace. Thus, the fields around Herbig Ae/Be 
stars represent ideal targets to study the transition from 
aggregates to dense clusters and to empirically probe the 
onset of cluster formation. 



Multiwalength studies of the environments of several 
Herbig Ae/Be stars in the optical, near infrared and mil- 
limet er (L kHa 101 

sony |1994|; BD+40°4124: Hillenbrand et al 
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et al. |1995| ; MacH12: Plazy & Menard |997| have shown 
that the young clusters are still partially embedded in the 
parent molecular clouds and that near infrared observa- 
tions, especially at K-band (2.2 ^m), are best suited to 
detect the less massive companions to the Herbig Ae/Be 
star itself. Observations at near infrared wavelengths of 
a consistent number (16) of fields around Herbig Ae/Be 
stars have been presented by Li et al. (1994). Their study 
is focused at the detection of bright companions and/or 
diffuse emission which may affect large beam photometry 
of the Herbig stars, and, for this reason, it is limited to a 
small area around each star, much smaller than the extent 
of the known clusters (e.g. LkHa 101, BD+40°4124). 

The first statistical study aimed at investigating the 
properties of star formation around Herbig Ae/Be stars 



using K-band images has been carried out by Hillenbrand (1995) 
In spite of the relatively small number of observed fields 
(17), Hillenbrand has found evidence for a correlation be- 
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tween the mass of the Herbig Ae/Be star and the surface 
density of K-band stars detected aro und it . These initial 
results were confirmed by Testi et al. ( 1997 ; hereafter Pa- 
per I), who obtained J, H and K images of moderately 
large fields around 19 Herbig Ae/Be stars. Using various 
methods to define the richness of the star cluster, in ad- 
dition to simple K-band star counts, in Paper I we have 
shown a clear dependence of the richness of the embedded 
clusters with the spectral type of the Herbig Ae/Be star. 
Moreover, our data seem to indicate that the clustered 
mode of star formation appears at a detectable level only 
for stars of spectral type B5-B7 or earlier. Whether the 
transition from isolated stars to clusters has a threshold 
or is a smooth function of stellar type (or, better, mass) 
is still unclear because of the small statistical significance 
of the samples studied so far. 

In order to complete our systematic study of Herbig 
stars, we have observed in the near infrared a second sam- 
ple of 26 fields selected from the compilations of Finken- 



zeller & Mundt Ql984| ) and The et al. ( |1994| ). Among the 
observed stars, we included Z CMa whose spectral type is 
F5 and whose membership to the Herbig group has been 
disputed. However, we will not consider it in the analy- 
sis of the clustering properties. Together with the 19 stars 
analysed in Paper I, our final sample consists of 44 objects 
covering almost uniformly the whole spectral range from 
09 to A7 stars. 

We have included in our sample 33 out of the 39 stars 
(85%) with decli nation greater than —11.5° listed in Finken- 
zeller & Mundt ( 1984 ) and replaced the remaining six ob- 
jects with 11 stars of similar sp ectra l type taken from the 
updated catalog of The et al. ( 1994 , Table I of members 
and candidate members). Thus, we are confident that the 
final sample gives a fairly complete representation of the 
whole class of Herbig Ae/Be stars and that the inferences 
discussed in this paper should have a solid physical and 
statistical meaning. 



In a companion paper (Testi et al. |1998| ; hereafter Pa- 
per II), we have collected the observational material (im- 
ages, colour-colour diagrams and stellar density profiles) 
of each of the combined sample of 44 fields. In this paper, 
we present the results of the analysis of this large body of 
observations aimed at the detection, characterization and 
comparison of the small star clusters around intermediate- 
mass PMS stars. 



2. Results 

2.1. Richness Indicators 

In our previous studies we have identified the most suit- 
able indicators of the richness of embedded stellar cluster 
that greatly reduce the problem of background/foreground 
contamination. The two quantities of interest are the num- 
ber of stars in the K-band image within 0.21 pc from 
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Fig. 1. Distribution of the cluster radius of 21 stellar 
groups detected by computing the source surface density 
profiles obtained in Paper II. 



the Herbig star with an absolute K magnitude Mk < 5.2 
(A/k), and the integral over distance of the source surface 
density profile subtracted by the average source density 
measured at the edge of each field, (Ic)- The choice of 
a radius of 0.2 pc for computing Ak is suggested by the 
typical value of the cluster size determined in Paper II, 
as the distance from the Herbig star where the sources 
surface density profile reaches the background level. As 
illustrated in Fig. |], the distribution of the radius of the 
stellar density enhancement shows a clear peak around 
r ~0.2 pc. This result is in good agreement with that 



found in various young stellar clusters (Hillenbrand 1995 



Carpenter et al. 1997). It is remarkable that stellar groups 
with a few to several hundred members share similar sizes, 
corresponding to the typical dimensions of dense cores in 
molecular clouds. 

The values of A/k and Ic computed for the 44 stars 
of our complete sample are given in Columns 3 and 4 of 
Table [l], together with the name and spectral type of the 
Herbig Ae/Be star at the center of each field. The uncer- 
tainty on Ic has been calculated by propagating the error 
in the determination of the background stellar density at 
large distances from the Herbig star. As in Paper I, we 
discuss the dependence of the richness indicators on the 
spectral type of the Herbig star rather than on its mass, 
even though the latter is the relevant physical parame- 
ter. The main reasons for this choice are the proximity of 
the Herbig star to the main sequence (and hence a close 
relation between spectral type and mass is "almost" well 
defined) and the much larger uncertainties in the mass 
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Table 1. Values of the richness indicators. 



Star 


Type 


Age (Myr) 


Ak 


Ic 


Ak 


A/g 


Til 

Ic 


Ic 


Log(pAA K ) 


Log(pi c ) 


V645 Cyg 


07 


— 


>5 


29.5 ± 2 


>5 


>5 


> 29.5 ± 2 


> 29.5 ± 2 


2.1 


2.9 


MWC 297 


09 


— 


37 


20.4 ± 1 


24 


23 


14.3 ± 1 


13.3 ± 1 


3.0 


2.8 


MWC 137 


B0 


— 


>59 


76.0 ± 9 


57 


55 


70.1 ± 5 


64.4 ± 4 


3.2 


3.4 


R Mon 


BO 


— 





-12.8 ± 3 








-5.1 ±3 


-5.1 ±3 


0.0 


0.0 


BHJ 71 


B0 


— 


4 


4.0 ±3 


1 


1 


-0.4 ± 1 


-0.3 ± 1 


2.0 


2.1 


MWC 1080 


BO 


— 


>9 


31.0 ±3 


>9 


>9 


> 31.0 ±3 


> 31.0 ±3 


2.4 


3.0 


AS 310 


B0 


— 


>37 


70.0 ± 17 


>34 


>34 


> 70.2 ± 17 


> 70.2 ± 17 


3.0 


3.3 


RNO 6 


Bl 


— 


>11 


11.0 ± 1 


>11 


>11 


> 11.0 ± 1 


> 11.0 ± 1 


2.5 


2.5 


HD 52721 


B2 


— 


10 


20.5 ±4 


5 


5 


11.1 ± 8 


11.1 ± 8 


2.4 


2.8 


BD+65°1637 


B2 


— 


29 


75.0 ± 5 


24 


28 


58.4 ± 4 


63.4 ± 4 


2.9 


3.3 


HD 216629 


B2 


— 


29 


34.0 ± 6 


23 


22 


27.0 ± 6 


25.3 ± 7 


2.9 


3.0 


BD+40°4124 


B2 


— 


19 


11.0 ±3 


16 


16 


12.6 ± 2 


12.9 ± 2 


2.7 


2.5 


HD 37490 


B3 


— 


9 


9.9 ±3 


6 


6 


3.4 ± 1 


3.8 ± 1 


2.4 


2.5 


HD 200775 


B3 


— 


8 


1.9 ± 1 


7 


7 


1.0 ±1 


1.0 ±1 


2.3 


1.8 


MWC 300 


Be 


— 


>2 


21.0 ± 8 


>2 


>2 


> 21.9 ±8 


> 21.9 ±8 


1.7 


2.8 


RNO IB 


Be 


— 


12 


9.7 ± 1 


12 


12 


8.9 ± 1 


9.0 ±8 


2.5 


2.5 


HD 259431 


B5 


0.05 


2 


0.9 ±2 








-1.7 ±3 


-0.2 ±2 


1.7 


1.4 


XY Per 


B6 


2.0 


3 


11.3 ± 3 


2 


5 


-0.5 ± 1 


-0.1 ± 1 


1.9 


2.5 


LkHa 25 


B7 


0.2 


11 


14.5 ± 5 


8 


5 


10.3 ± 3 


9.0 ±3 


2.5 


2.6 


HD 250550 


B7 


0.5 


4 


2.2 ± 2 


2 


2 


0.8 ± 1 


0.8 ± 1 


2.0 


1.8 


LkHa 215 


B7 


0.1 


7 


3.9 ± 1 


6 


4 


5.4 ± 1 


3.6 ± 1 


2.3 


2.1 


LkHa 257 


B8 


2.0 


15 


5.5 ±6 


14 


16 


8.3 ±3 


9.5 ±4 


2.6 


2.2 


BD+61°154 


B8 


0.1 


8 


-1.4 ±3 


4 


4 


2.0 ± 1 


2.4 ± 1 


2.3 


0.0 


VY Mon 


B8 


0.05 


25 


23.2 ± 5 


18 


16 


14.9 ± 2 


12.9 ± 2 


2.8 


2.8 


VV Ser 


B9 


1.0 


24 


16.9 ± 5 


20 


23 


3.1 ± 2 


3.0 ± 2 


2.8 


2.7 


V380 Ori 


B9 


1.0 


3 


-2.0 ±2 


3 


3 


-1.5 ± 3 


-1.5 ± 3 


1.9 


0.0 


V1012 Ori 


B9 


— 


4 


1.9 ± 2 


4 


4 


1.0 ± 1 


1.0 ± 1 


2.0 


1.8 


LkHa 218 


B9 


1.0 


8 


2.0 ±5 


5 


7 


1.7 ± 1 


2.0 ± 1 


2.3 


1.8 


AB Aur 


AO 


1.5 


>3 


3.0 ±6 


>1 


>3 


-1.0 ± 2 


-2.3 ± 2 


1.9 


2.0 


VX Cas 


AO 


1.0 


13 


4.5 ±4 


34 


35 


3.9 ±4 


4.4 ±4 


2.5 


2.1 


HD 245185 


A2 


7.0 


10 


4.5 ±5 


19 


19 


5.0 ± 1 


5.0 ± 1 


2.4 


2.1 


MWC 480 


A2 


7.0 


>3 


5.0 ±6 


>9 


>9 


-2.4 ±2 


-2.4 ±2 


1.9 


2.2 


UX Ori 


A2 


1.0 





-0.3 ± 1 








-0.1 ± 1 


-0.1 ± 1 


0.0 


0.0 


T Ori 


A3 


1.0 


5 


1.0 ±2 


4 


7 


0.4 ± 1 


1.9 ± 2 


2.1 


1.5 


IP Per 


A3 


7.0 


3 


5.3 ±4 


6 


6 


-0.7 ±3 


-0.7 ±3 


1.9 


2.2 


LkHa 208 


A3 


1.0 


4 


2.2 ±5 


3 


2 


1.2 ± 1 


2.0 ± 1 


2.0 


1.8 


MWC 758 


A3 


6.0 


>2 


3.4 ± 1 


>3 


>3 


-0.3 ± 1 


-0.3 ± 1 


1.7 


2.0 


RR Tau 


A4 


0.1 


7 


0.8 ±6 


2 


2 


-0.4 ± 4 


-3.2 ± 3 


2.3 


1.4 


HK Ori 


A4 


7.5 


7 


2.2 ± 1 


11 


11 


2.2 ± 1 


2.2 ± 1 


2.3 


1.8 


Mac H12 


A5 




15 


5.1 ± 1 


21 


21 


5.9 ±1 


6.6 ±1 


2.6 


2.2 


LkHa 198 


A5 


10.0 


6 


-10.6 ± 11 


14 


14 


-9.4 ± 10 


-8.4 ±9 


2.2 


0.0 


Elias 1 


A6 




>2 


2.0 ±3 


>2 


>2 


1.0 ±1 


1.0 ±1 


1.7 


1.8 


BF Ori 


A7 


3.0 


4 


1.1 ± 1 


3 


3 


1.1 ± 1 


1.1 ± 1 


2.0 


1.5 


LkHa 233 


A7 


4.0 


2 


1.0 ± 1 


4 


5 


0.6 ± 1 


1.5 ±1 


1.7 


1.5 



estimates based on the location of the stars in the HR di- 
agram and the use of evolutionary tracks. On the other 
hand, spectral types of young stars sample are usually de- 
termined with an uncertainty of one or two subclasses, 
with the exception of a few cases (see the discussion in 
Paper II). In particular, RNO IB and MWC 300 do not 
have a reliable classification and in the literature they are 
generically referred to as Be stars. For plotting purposes, 
we have arbitrarily assigned a B5 classification to both of 
them. These stars have both values of Ic and A/k typical 



of B stars, and the assigned spectral type of B5 does not 
affect our results. 

The variation of A/k and Ic as a function of the spec- 
tral type of the Herbig star is shown in Figures ^ and |[ 
respectively. Both figures confirm the initial results found 
in Paper I and by Hillenbrand 1995 that higher mass stars 



tend to be surrounded by richer clusters. As expected, the 
evidence for clusters is more pronounced using Ic instead 
of A/k , but the fact that the same trend is seen in both 
indicators gives strength to its real significance. 
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B5 AO 

Sp. Type 



FO 




B5 AO 

Sp. Type 



FO 



Fig. 2. A/k as a function of the spectral type of the Herbig Fig. 3. Ic versus spectral type for the whole sample. As 
Ae/Be star. The two stars with an uncertain spectral type in Fig. ||, the two stars with an uncertain spectral type Be 



Be have been plotted as B5. 



have been plotted as B5. 



In addition to the evidence for a variation of Ic with 
spectral type, the results of Fig. || suggest the existence 
of three regimes for the distribution of stars around Her- 
big stars. In the first one, characterized by Ic ~ 40, the 
Herbig stars are definitely associated with rich clusters; in 
the intermediate regime with 10 $ Ic~40 a small cluster 
may be present; in the third case where Ic ~ 10 only small 
aggregates or background stars in the field are found, a 
situation similar to that observed in low-mass star form- 
ing regions. Three stars of our sample (MWC 137, AS 310 
and BD+65°1637) belong to the first group. The sizes of 
the clusters derived for these stars in Paper II are on the 
large side of the distribution (~ 0.4 pc), however, at least 
in two cases (AS 310 and BD+65°1637) the Herbig star is 
not at the center of the cluster, and this results in an over- 
estimate of the cluster radius, determined from the radial 
density profile. Among stars of early spectral types, we 
find a large spread of values of Ic- In some cases, the low 
I c are dubious. For example, V645 Cyg and MWC 1080, 
have reletively low values of Ic typical of the intermedi- 
ate regime. However, both stars are embedded in bright 
nebulosities which may affect the star count resulting in a 
severe underestimate of the actual stellar density. In other 
cases, however, the low values of Ic may be real. 

The presence of groups or small clusters is secure around 
most of the stars belonging to the intermediate regime. A 
clear central density peak is observed in the density pro- 
files of almost all stars of this group, as shown in Paper II. 
Typically, the Herbig star is located at the center of the 
stellar group, which generally has a round shape. A pos- 
sible exception is the VY Mon field, where the stellar ag- 



gregate has an elongated structure with an aspect ratio 
of approximately 4:1 (the Herbig star is at the center). 
HD 52721 is not associated with molecular gas (Fuente 
et al. 1998), it is a rather old star, and the large cluster 
radius is probably the result of dynamical relaxation. 

Finally, the large majority (~ 65%) of the Herbig stars 
of the sample have very low values of Ic showing no 
enhancements above the background stellar density. All 
the fields around stars with spectral type later than B9 
belong to this group. However, although most of these 
stars have spectral types later than B7-B8, there are ex- 
treme cases of early-type stars that deserve some discus- 
sion. First, the negative value of Ic found for R Mon 
(B0) is probably due to localized extinction around the 
star, which is probably on the observer side of a molecu- 
lar clump. Also, the stars BHJ 71 (B0), HD 200775 (B3) 
and HD 259431 (B5) have anomalously low values of Ic- 
The first one is a poorly studied star, and we have not 
been able to find any information on the molecular gas 
in the literature (see Paper II). HD 200775 is in many 
respect very similar to BD 65° 1637: same spectral type, 
same molecular gas morphology and similar evolutionary 
status (Fuente et al. 1998 ). The absence of a cluster of 
young stars around HD 200775, as opposed to the rich 
cluster around BD 65° 1637, may be the result of dynam- 
ical dissipation (see below) or reflect a difference in the 
environment. HD 259431 appears to be isolated in spite 
of the large amount of molecular material around it (Hil- 
lenbrand 1995). Among late- type stars, LkHa 198 stands 
out for its negative value of Ic due to the bright nebu- 
losity and localized extinction. Interesting cases with an 
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Fig. 4. A/jj an d A/k as a function of the spectral type of the 
Herbig Ae/Be star. Open circles denote fields for which we 
could not compute the age using the Herbig Ae/Be star 
parameters (see Paper II and sect. 2.2 of the main text). 
As in Fig. ^, the two stars with an uncertain spectral type 
Be have been plotted as B5. 

indication of an extended population of embedded sources 
but not spatially concentrated around the Herbig star are 
VX Cas (AO), HD 245185 (A2) and Mac H12 (A5). 

2.2. Mass sensitivity correction 

Since young stars change their bolometric luminosity and 
effective temperature during PMS evolution, the mass of 
the smallest star detectable (Mi) in our K-band images 
is a function of age (t c ), of the absolute completeness K- 
magnitude (M^) and of the extinction along the line of 
sight (A K ). 

In Paper II, using the PMS evolutionary tracks for in- 
termediate mass stars of Palla & Stahler (1993) and the 
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Fig. 5. 1^ and Iq versus spectral type for the whole sam- 
ple. Open circles as in Fig. |[ As in Fig. |[ the two stars 
with an uncertain spectral type Be have been plotted as 
B5. 



compilation of stellar parameters from the literature, we 
have computed in an homogeneous way the ages of most of 
the stars later than B5. Then, using the PMS evolut ionar y 
tracks for low mass stars of D'Antona & Mazzitelli ( 1994 ), 
we have translated M^- into Mi for two values of the ex- 
tinction (Ak = and 2) assuming that all the stars in 
each group are coheval with the Herbig Ae/Be star. 

There are some stars (mostly early Be types) for which 
it is not possible to derive ages from PMS tracks. For these 
stars, we have estimated the age in the following way. From 
the results of Paper II we see that the Herbig Ae systems 
tend to be 1 to 10 Myr old (the mean value is ~ 4.6 Myr), 
while Be systems tend to be younger than 1 Myr (with a 
mean age of ~ 0.7 Myr). Thus we decided to adopt as age 
t c = 0.5 Myr for the 17 Be systems and t c = 5 Myr for the 
2 Ae systems without an age determination in Paper II. 
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If Ak ~0 mag, all the fields with a good age estimate 
from Paper II are complete to less than 0.1 Mq, with 
the only exception of the two fields around HK Ori and 
LkHod98. In this case, A/"k and I c sample the totality of 
the stellar population in practically all fields. If Ak ~ 
2 mag, the minimum mass is generally larger, ranging 
from < 0.1 to 0.49 M Q in HK Ori. For each field we have 
computed the absolute magnitude corresponding to these 
masses and the corresponding extinctions at K-band, then 
we have calculated the two richness indicators discussed in 
Sect. ^ considering, in each field, only the stars with mag- 
nitude lower than the computed one. In this way, under 
the assumption that the mean extinction is approximately 
the same in all fields (and for all stars in each field), we 
have calculated mass limited richness indicators. The val- 
ues of A/k and I c for A K = and 2 (A/^,A/"|, I c and I c ) 
are reported in Table [|, Columns 6-9. Since some of the 
fields without an estimate of the age of the Herbig Ae/Be 
star have mass sensitivities lower than our limits (which 
are based on HK Ori), some of the values reported are 
lower limits (indicated by >). 

In Fig. | and | we show the behaviour of the "mass 
sensitivity corrected" indicators. Fig. [| clearly shows that 
trend detected in Fig. || is completely cancelled by the 
bias introduced by the contamination from field stars. In 
fact, since the late type stars are older than early type the 
absolute magnitude corresponding to the same mass limit 
is higher for the Ae systems than for the Be systems, and 
since this discrepancy is not compensated by the difference 
in distance, we are effectively counting more field objects 
in Ae fields than in Be fields. On the contrary, Fig. || 
clearly shows the same trend as Fig. |^ supporting our 
conclusion of a dependence of Ic on the spectral type of 
the Herbig Ae/Be star. Therefore in the following we will 
use Ic which has a higher signal to noise ratio than I c and 
T 2 

3. Discussion 

In Paper I we suggested the presence of a threshold for 
the appearence of clusters around Herbig stars of spec- 
tral type earlier than ~B7. The results discussed in the 
previous section confirm the property that Be systems are 
associated with more conspicuous groups than the Ae sys- 
tems. Although the identification of a critical spectral type 
(or mass) is no longer evident in the complete sample, it 
is worth noting that no stars later than B9 show large 
groups, i.e. their Ic is < 10. 

This general result raises some important issues: (1) is 
this difference between Be and Ae system related to a dif- 
ferent location in the galaxy? (2) since Ae systems are gen- 
erally older than Be systems, could the disappearence of 
the cluster be an evolutionary effect? (3) what is the typ- 
ical environment in which intermediate mass stars form? 
(4) does the standard accretion scenario for the forma- 



tion of Herbig stars need to be replaced by formation in 
clusters? 

The first question is easy to answer. We do not find 
evidence for a dependence of the clustering properties with 
the galactic position of the target star, shown in Paper II. 
We must caution, however, that most of the observed stars 
lie in the outer regions of the Galaxy, due to the selection 
effect introduced by our observing sites in the northern 
hemisphere. Although we do not expect an opposite result 
toward the inner Galaxy, we must await the results of a 
similar study on a southern sample before generalizing our 
conclusion. 



3.1. Age effect 

Following the age estimates given in Paper II, our sample 
of Herbig Ae/Be stars covers the range of ages from ~ 0.05 
to 10 Myr. As noted above, Ae stars tend to be older 
than Be stars. Thus, it is possible that the variation of 
the richness indicators from Be to Ae systems is caused 
by some evolutionary effect. It is possible that a stellar 
group composed by few tens or hundreds of objects could 
be dispersed on a timescale of a few million years and 
become undetectable (i.e. confused with the field stars) in 
older systems. 

In order to explore this possibility more quantitatively, 
we show in Fig. |^ the run of Ic as a function of the age 
for the 25 fields with an age determination of the Herbig 
star (see Table 2 of paper II) . The lack of stars with high 
values of Ic for ages greater than 2 Myr is related to the 
fact that our sample does not contain Be stars of that 
age. However, it is worth noting that for Be systems (filled 
circles in Fig. ^) younger than 2 Myr we find high and low 
values of Ic at any age. 

From the figure we see that it does not seem to exist 
a correlation between the age of the Herbig star and the 
presence of a cluster. The correlation is only with the spec- 
tral type (or mass). Stars earlier than B7 have clusters or 
not independently of their age; stars later than B7 do not 
have clusters and are on average older than the rest. If we 
say for convenience that a cluster is detected for Ic Z 10, 
the oldest star surrounded by a cluster is XY Per with an 
age of 2 Myr. Since stars of age greater than ~4 Myr do 
not have clusters, there are two possibilities to account for 
such property: these stars did form in relative isolation or 
the attending cluster has by now disappeared. The former 
explanation is obvious and does not require any further 
comment. 

The second alternative is theoretically more interest- 
ing and we have performed several N-body simulations to 
verify under which conditions the dynamical evolution of 
a small cluster can evolve so rapidly to dissipate the ma- 
jority of its members in a time scale less than ^4 Myr. We 
have considered models containing N =100, 150 and 200 
stars starting from different initial conditions (virialized 
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tion to the upper limit of the mass distribution. Only if 
the most massive star exceeds 5 Mq, the evolution is fast 
enough that the system loses many members in several 
crossing times. Otherwise, it is impossible to modify the 
cluster composition in a few million years. Now, since a 
5 M Q corresponds to a B6 ZAMS star, it is clear that 
Herbig stars of later types which were formed in relatively 
small clusters would have retained the original population 
of companions at the observed age. Therefore, we conclude 
that the fact that we do not see evidence for the presence 
of clusters around stars of type later than about B7-B8 is 
an imprint of the stellar formation mode rather than the 
consequence of the dynamical evolution of a presently dis- 
persed cluster. On the other hand, the absence of clusters 
around some of the early type B stars could result from 
the dynamical evolution, especially in case of a relatively 
small initial population. 



3.2. The transition from loose to dense groups 



Fig. 6. Ic versus age. Only the fields with age determi- 
nation from Paper II are presented. The arrow represents 
the field around LkHa 198 (Ic = -10.6). filled circles: Be 
systems; open circles: Ae systems. 



and non virialized clusters) and with varying frequency 
distributions of stellar masses (Salpeter, Scalo, Kroupa 
and uniform mass). The half-mass radius was varied be- 
tween 0.4 and 0.7 pc and the typical velocity dispersion 
was 2 km s . Thus, the clusters are initially richer and 
larger than those actually observed around Herbig stars 
characterized by Ic ~30-40 and a radius of ~0.2 pc. We 
have also varied the upper limit of the most massive ob- 
ject from 20 Mq, corresponding to a 09 star, to 2 Mq, 
corresponding to an A5 star. It is well known that the 
dynamical evolution of cluster models is very dependent 
on the choice of the IMF, especially in what concerns the 
effects of mass segregation at early times (e.g. de la Fuente 
Marcos 1995, 1997). We have not considered the presence 
of primordial binaries in the initial population. 

For most runs, we have found that it is very difficult 
to loose a significant fraction of stars in the initial few 
crossing times, as required by the age constraints of the 
Herbig stars. The only cases where the stellar population 
decreases to less than half the initial value in about ten 
crossing times are those characterized by the lowest num- 
ber of stars (N =100 in our models) distributed in mass 
according to a Salpeter IMF. The first point is a criti- 
cal one, since the evaporation time increases very rapidly 
with the number of stars. Even models with N = 150 do 
not evolve sufficiently rapidly. As for the IMF, the result 
is acceptable even though we know that a single power- 
law approximation must break down at subsolar values. 
More important, however, is the sensitivity of the evolu- 



The values of A/k and Ic can be transformed into star 
(or effective star) volume densities using the average clus- 
ter radius of 0.2 pc and assuming spherical clusters. The 
logarithm of the stellar densities in sources per cubic par- 
sec are reported in the last two columns of Table [l| This 
estimates of the stellar densities are affected by large un- 
certainties (as high as 30%, Hillenbrand 1995), due to the 
unknown depth of the clusters along the line of sight and 
the neglection of concentration effects, but it is a useful ap- 
proximation to be used for comparing the various regions. 
The resulting distributions shown in Fig. ^ confirm the 
apparent trend discussed in Sect. 2 and put on a firmer 
ground the suggestion of Hillenbrand (1995) of a physi- 
cal relationship between cluster density and the maximum 
mass of the Herbig star. For comparison, we also show in 
Fig. |^ the stellar density range of the T Tauri aggregates 
found in Taurus- Auriga by Gomez et al. ( 1993 ) (shown as 
the solid line at G/K, not to scale), and the same qu antity 
for the Trap ezium cluster (McCaughrean & Stauffer 1994 ; 
Hillenbrand 1997 ). The conclusion is evident: intermedi- 
ate mass stars mark the transition from low density ag- 
gregates of < 10 stars per cubic parsec of T Tauri stars to 
dense clusters of > 10 3 stars per cubic parsec associated 
with early-type stars. 



3.3. Implications for star formation 

The conventional theory of protostellar infall successfully 
accounts for the formation and early stellar evolution of 
low- and intermediate-mass stars up to about 10 Mq. 
The location of the stellar birthline and the distribution 
of the observed T Tauri and Herbig Ae/Be stars in the 
HR diagram are in excellent agreement with the theoreti- 
cal predictions (Palla & Stahler 1990, 1993). Protostars 
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Fig. 7. Stellar volume densities derived from A/k (fe/t) and from Ic (right) versus spectral type of the central star. 
Stars with Ic < have been excluded. The heavy vertical line at 06 represents the range of stellar densities found in 
the Trapezium cluster, whereas that at G/K (not to scale) represents the densities of stellar groups in Taurus- Auriga. 



more massive than about 10 M Q burn hydrogen while 
still in the accretion phase and therefore join the main 
sequence, implying that stars of even higher mass have no 
contracting PMS phase. Observations of clusters such as 
NGC 6611 and the dense regions of the Trapezium cluster 
support this finding (e.g. Hillenbrand 1993, 1998). How- 
ever, the accretion scenario of isolated protostars fails to 
explain the existence of massive stars. Radiation pressure 
by photons produced at the stellar and disk surfaces on 
the infalling gas begins to become significant at about the 
critical mass of 10-15 M (e.g. Larson & Starrfield 1974; 
Yorke & Krugel 1977). This limit can be increased by con- 
sidering variations in the dust properties (abundances and 
sizes) or in the mass accretion rate (Wolfire & Cassinelli 
1987). But in either case the required conditions are so 
extreme (dust depletion by at least one order of magni- 
tude and accretion rates of at least 10 -3 M yr _1 ) that 
cannot reasonably apply in all circumstances. Any depar- 
ture of the infall from spherical symmetry can also help to 
shift the critical mass to very massive objects (e.g. Nakano 
1989; Jijina & Adams 1996). 

The fact that, as we have seen, stars form in groups 
and clusters rather than in isolation may offer an alternate 
mechanism to circumvent the problem of the high lumi- 
nosity of massive stars and the negative feedback on the 
accreting gas. Accretion processes in clusters has long been 
advocated to explain the shape of the initial mass function 
(Zinnecker 1982; Larson 1992; Price and Podsiadlowski 
1995), and to account for the location of massive stars 
in the cluster centers (Bonncll et al. 1997). Very recently, 
Bonnell et al. (1998) have presented model calculations of 
the dynamical evolution of rich, young clusters that show 



that accretion-induced collisions of low- and intermediate- 
mass stars (formed in the standard accretion mode) can 
result in the formation of more massive objects in time 
scales less than 10 6 years. However, the process requires 
a critical stellar/protostellar density so that collisions and 
merging are frequent. For the formation of a star of mass 
of ~50 M , an initial density of >10 4 stars pc -3 is re- 
quired. This is the typical density observed in the central 
regions of large clusters such as the Orion Nebula clus- 
ter where stars so massive are indeed present, but it is 
hardly the case for the stars of our sample. As shown in 
Fig. 0, the highest stellar volume densities are limited to 
about 2xl0 3 stars pc~ 3 , almost an order of magnitude 
lower than the minimum value for the efficient collisional 
build-up of massive stars. These densities are found for 
stars of spectral type B0-B2, i.e. with ZAMS masses of 
about 20 to 10 M . Thus, these stars are at the border- 
line between the conditions of isolated and collisional ac- 
cretion mechanisms. It would be extremely interesting to 
extend the Bonnell et al. (1998) models to less extreme 
conditions on the stellar density in order to explore the 
dynamical evolution of clusters resembling those observed 
around Herbig Be stars. Based on the observational evi- 
dence, we may conclude that in most cases, the formation 
of Herbig stars can be understood in terms of the conven- 
tional accretion scenario, whereas for the most massive 
members of this group dynamical interactions in the clus- 
ter core and residual gas accretion can be of importance 
in determining the final mass, even though the observed 
stellar density is never extremely high. 
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4. Conclusions 

Using the techniques described in Paper I, we have ana- 
lyzed the near-infrared observations of the fields surround- 
ing Herbig Ae/Be stars presented in Paper II with the 
goal of identifying and characterizing the presence of (par- 
tially) embedded clusters formed around the intermediate- 
mass stars. We have examined 44 fields around stars rang- 
ing in spectral types from A7 to 09. The main results of 
these studies can be summarized as follows. 

We have confirmed and extended the correlation be- 
tween the spectral type (or mass) of the Herbig stars and 
the number of nearby, lower mass objects, first noted by 
Hillenbrand (1995). Rich clusters with densities up to 10 3 
pc~ 3 only appear around stars earlier than B5, corre- 
sponding to a mass of about 6 M Q . Conversely, A- type 
stars are never accompanied by conspicuous groups and 
the typical density is less than 50 pc~ 3 . However, the 
transition between formation in loose groups and in clus- 
ters does not occur sharply around spectral type B7 as 
we suggested in Paper I. The appearence of denser stellar 
groups is quite smooth moving from Ae to Be systems, 
thus suggesting that intermediate-mass stars fill naturally 
the gap between the low-density, low-mass aggregates and 
the high-density, high-mass clusters. 

Using a richness indicator, Ic, based on stellar sur- 
face density profiles, we have identified three regimes for 
the distribution of stars around Herbig Ae/Be stars: rich 
clusters characterized by values of Ic i 40; small clusters 
or aggregates with 10<Ic~40, and small aggregates or 
background stars for lower values of Ic ■ Herbig stars with 
rich clusters are rare: only three stars of the sample defi- 
nitely belong to the first group, whereas in other cases we 
have indication that the star counts may severely under- 
estimate the actual stellar density. The typical cluster size 
is ^0.2 pc irrespective of the richness of the clusters. This 
value is remarkably close to the dimension of dense cores 
in molecular clouds. 

The majority of Herbig stars (65%) are found in the 
last regime, showing no enhancement above the background 
stellar density. All the stars with spectral types later than 
B9 belong to this regime and are generally older than Be 
stars. This result is not affected by different sensitivities 
to the lowest masses in different fields. From dynamical 
considerations, we conclude that the absence of clusters 
around late-Be and Ae stars is an imprint of stellar for- 
mation in relatively isolated cores (as in the case of T 
Tauri stars), and not the result of rapid cluster evolution 
and dispersion. 

The fact that the most massive Herbig stars have clus- 
ters strongly supports the notion that their formation has 
been influenced by dynamical interactions with lower mass 
stars and/or protostellar cores. It is no coincidence that 
the onset of clusters manifests at a mass of about 8-10 M 
where the conventional accretion scenario for isolated pro- 
tostars faces severe theoretical problems. Future studies of 



the formation and evolution of these stars should take into 
account the observational evidence for high stellar density 
environments. 
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